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Abstract 

Localization of flux lines to splayed columnar pins is studied. A sine-Gordon 
type renormalization group study reveals the existence of a Splay glass phase 
and yields an analytic form for the transition temperature into the glass phase. 
As an independent test, the I—V characteristics are determined via a Molecu- 
lar Dynamics code. The glass transition temperature supports the RG results 

convincingly. The full phase diagram of the model is constructed. 
PACS numbers: 74.20.Mn, 74.60.Ge, 64.60.Ak, 74.40. +k 
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Zero resistance in Type II superconductors is maintained in high magnetic fields by 
pinning the flux lines to avoid dissipation. The central issue of technological interest is how 
to maximize the critical current J c which depins the vortices. Atomic vacancies or light 
particle irradiation typically modelled by random point disorder collectively pin the flux 
lines into an entangled arrangement leading to the formation of the Vortex glass phase [1]. 
On the other hand, columnar pins produced for example by heavy ion irradiation, twinning, 
or screw dislocations give rise to the Bose glass phase [2,3]. It was argued that J c is strongly 
enhanced by columnar pinning over the traditional point defects [2]. The two types of 
glasses are characterized by distinct localization mechanisms, manifested in different I-V 
curves. Recently the phase diagram of the model with the simultaneous presence of both 
types of disorder has been constructed [4]. The main result was that the nature of the 
localization is dominantly of the Bose glass type with a reduced transition temperature [5]. 

In the quest of maximizing J c , the idea of tilting the columnar defects was proposed [6]. 
In the case of parallel columns, once a segment of the vortex reaches an adjacent column 
either via thermal activation or a quantum mechanical process, the remaining part of the 
vortex can follow at no additional energy cost. In the Splay glass phase, during a vortex 
hop an ever increasing segment of the line will be forced into an energetically unfavorable 
region. This provides an effective confining mechanism and leads to a dramatically enhanced 
value of J c . Subsequent experiments have indeed confirmed this expectation [7]. In spite 
of intense interest, the phase diagram of the model has been determined only qualitatively. 
The main result of the present work is the quantitative determination of the phase diagram. 
First a renormalization group (RG) analysis is used to determine the dependence of T c on 
the degree of columnar tilting and flux line rigidity. This RG treatment is then shown to 
be in excellent agreement with molecular dynamics simulations. From a theoretical point of 
view the problem is also interesting because of its relevance to Quantum Phase transitions 
occuring in boson systems [3,8]. 

Our analysis begins with the Hamiltonian for N interacting flux lines confined to a plane, 



2 



N i 1 r 7 / \ i 2 



dr t (z) 



dz 



+ jE%( 2 )-^)] + w),4 (i; 



Here the scalar displacement field r 4 -(z) describes the position of the i-th flux line in a 
direction perpendicular to the applied field (z), vq denotes the strength of the line-line 
interaction, V denotes the disorder potential, and ei the elastic energy cost per unit length 
for transverse wandering. We choose to work with this restricted model since it is more 
amenable to analytic and numerical treatment than the corresponding 3-D model. It has 
been recently suggested via scaling arguments that the fluctuations of the flux line along 
the direction of the applied current are irrelevant at the transition and thus a planar model 
may be immediately applicable to 3-D systems [9]. 

The correlator A of the disorder potential V is defined via 

{V(x, z)V{x', z')) = A(x -x',z- z'). (2) 

For the case of uncorrelated point disorder, A p (x,z) ~ 6(x)6(z). Columnar disorder is 
correlated along the ion trajectories. When aligned parallel to the applied magnetic field 
A c (x,z) ~ 6{x). These two cases have been analyzed separately in Refs. [1-3], while the 
competition between the two types of disorder was investigated in Refs. [4,5]. In the present 
work, the disorder correlator represents tilted columnar lines with a Gaussian distribution 
of tilting angles of width vjj. The Fourier transform of the correlator is 

A s (q X} q z ) = A r— ^ — - exp[-(q z A x ) 2 /2v 2 D (q x A z ) 2 ] 
V2,Trv D \q x \ 

= A h s (q x ,q z ), (3) 

as was also used in Ref. [6]. Here A XjZ are the ultraviolet cut-offs and we have neglected the 
width of the vortex core. For the case of vd = 0, A recovers the untilted columnar form, 
while vd — > oo describes an isotropic distribution of tilting angles. 
The Hamiltonian can be rewritten in a sine- Gordon-like form as 
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Here u is the displacement field in a continuum description. We have applied a coarse- 
graining procedure and introduced replicas (indexed by a, f3 = f,---,n) to handle the 
disorder. The parameters K x = Vq/T and K z = ati/T are the bulk and tilt moduli, 
respectively, which describe the elasticity of the flux line array, and = A /T 2 . The 
bare coupling constant of the nonlinear term, representing the interactions, is given by 
g = A /(aT) 2 , with a the average line separation. As shown in Refs. [4,5], the transforma- 
tion u(x,z) — > u(x,z) + f x dxV (x , z) before replicating Eq. (4) allows us to absorb the 
forward scattering in the x direction oc d x u at the expense of generating a forward scattering 
in the z direction oc d z u. We thus obtain 
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This has the appealing feature that in the case of no tilting (vd = 0), the Bose glass 
transition is immediately recovered [2,3]. 

We have performed a renormalization-group analysis of the model Eq. (4) to first order 
in the interaction g. We find the following recursion relations upon rescaling by a factor e l 
and taking the n — > limit: 
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Here K~ x = 2tt(K x K z )~ 1 / 2 ~ T and the ratio of the tilt to bulk modulus, a = K z /K x = 
€ia/p describes the elastic anisotropy of the flux array, Ci^ are positive constants, and / is a 
function to be specified below. The bulk modulus K x is unrenormalized at first order (unlike 
the case of point or columnar disorder, our model does not possess statistical invariance to 



tilt [4]), while the renormalization of the tilt modulus K z only occurs for no tilting. This 
distinguishes the Bose glass (vd = 0) from the Splay glass, (vd ^ 0). In the former, the 
tilt modulus scales to infinity, while in the latter it remains finite both above and below 
the transition. This suggests that the Bose glass is indeed unstable to the splaying of the 
columns and is replaced by a new thermodynamic phase. In fact, in many respects the 
recursion relations are similar to those obtained for the case of point disorder. Hence it is 
also termed "correlated Vortex glass" [6]. 

The most important feature of the recursion relations is the reduction of the glass tran- 
sition temperature not only for increased tilt of the columns but also for increased elastic 
anisotropy of the flux line array. The Splay glass transition temperature is given by 
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The function / can be written in terms of a confluent hypergeometric function U [10], 

f(x) = U(2,l/2,x/2). 

It has the limiting behavior f(x — s- 0) = 4/3 + 0(y/x), and f(x — > oo) = 4/x 2 — 40/x 3 + 
0(1/ x 4 ). For the case vd = 0, the previous results [2,3] for the Bose glass are recovered, 
Kg X G = 3. For finite tilting the localization transition occurs at a lower temperature, in 
agreement with the arguments made in Ref. [6]. The phase boundary is shown in Fig. (1) 
as a function of the width vd of the tilting angle distribution and the anisotropy parameter 
a. For small tilting we find, 

KsklK£ G . = 1 - 1 ^« 5/2 4 + 0(v* D ). (11) 

The transition temperature falls off more slowly with tilting angle (~ v^) than suggested 
from the arguments presented in Ref. [6]. From an application point of view, it is worth 
noting that the reduction of the transition temperature can be counteracted by increasing 
the elasticity anisotropy, as represented in Eq. (9) by the parameter a. This feature is absent 
for the case of straight columnar pins. Lastly we note that for fixed a, the glass transition 



temperature saturates at a non-zero value for the limit of vd — > oo, i.e. the isotropic limit. 
In this limit the columnar pins become more and more point-like and thus are relatively 
ineffective at pinning, leading to a correspondingly lower transition temperature. 

To check the predictions of the RG, we have performed molecular dynamics simulations 
of Eq. (I) directly. The most direct way to identify the glass phase is based on the study of 
the response of the flux lines to an external current [11,12]. Assuming overdamped dynamics, 
the equation of motion for the lines is given by 

d t r = rid 2 x r - d x (V p + V v ) + F + f. (12) 

Here rj sets the energy scale for bending of a line (77 oc ei), V PjV are the pinning and vortex 
interaction potentials, F is the Lorentz force exerted by the applied current, and / is the 
Langevin noise taken to model the effects of temperature T via 

{f(x, z } t)f(x', z', t')) = 2T6(x - x')6(z - z')6(t - t'). (13) 

The potentials are chosen as follows: (1) The disorder pinning potential is modeled as the 
smooth function V p (x } z) = —V p i n [(x / R p i n — l)(x/ R p i n + l)] 2 for | x \< R p i n and otherwise, 
where x measures the distance from the center of the pin, given by x = vz } to the line 
and v is the tangent of the angle of the splayed columns which is chosen randomly from a 
uniform distribution of width vjj. The position of the columns are randomly placed in the 
x — z plane. (2) The vortex interaction potential is taken to be V v = V vor t e x^~ x ^ Rvortex for 
I x |< 6R vor tex and otherwise, where x is the vortex- vortex separation. All length scales are 
measured in terms of R vor tex which is given approximately by the screening length A ~ 5000 
angstroms, while energy scales are set by V vorte x = 2($ /47rA) 2 , with $ the flux quantum. 
We lastly chose to work with narrow pinning defects so that R p i n = 0.3R vorte x- 

The parameters for our simulations were chosen as follows: our system size was 
32 X 32 in units of R vor t eX i an d contained 20 lines and 32 columns. The parameters 
i) = 2.5, Vpi n = V vor tex = 1 were fixed for all runs. While non-unique, this choice of the 
parameters falls within the range of values pertinent to cuprate superconductors [13,8]. The 
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average velocity of the flux line array was measured for different temperatures and averaged 
over a varying number of disorder configurations so as to keep the statistical fluctuations 
below half a percent. Using l/V vorte x to set the time scale, typically less than 2,000 cycles 
were sufficient to reduce the velocity fluctuations to values less than f percent, thus ensuring 
proper equilibration of the system. A typical series of results are shown in Fig. (2) in a 
log-log plot. Here the tilting angle distribution width was Vd = 2. Two regions can be 
clearly seen at low temperatures - a linear response at large driving forces corresponding to 
free flux motion, and a non-linear response corresponding to thermally assisted flux motion. 
The crossover between the two becomes broader as the temperature is increased and the 
Ohmic region becomes larger. 

We define the Splay glass transition temperature as the temperature at which the entire 
curve obeyed an Ohmic response [14]. Our results for different tilting angles are shown in 
Fig. (3). Other methods for establishing the critical temperature, such as identifying the 
point of maximum slope of the I — V curves or monitoring when the fluctuations of the 
velocity due to the disorder drop below 1 percent, yield the same results to within the error 
bars. The solid line displays the transition temperature obtained from the RG method. The 
nonuniversal Bose glass transition temperature, which is 3 in the RG treatment, depends 
on the microscopic details of the model. In our numerical study Tb.g. was determined from 
simulations with no tilting of the columnar defects. With the above choice of parameters we 
obtained K^ G = 2.35 ±0.05. The solid line in Fig. (3) is subsequently obtained by inserting 
this value into Eq. (10). No adjustable parameters enter the expression. We have explored 
a larger region of the parameter space by changing the time step dt, rj, V vorte x /V P im and the 
line density, and found that these changes modified the overall transition temperature only 
through the combination defined by a as suggested by the functional form obtained from 
the RG treatment. The numerical data support the RG results remarkably well lending 
credence to the full phase diagram of the model given in Fig. 1. This is the central result 
of our paper. 

In conclusion, we have investigated the localization transition of flux lines in the presence 
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of tilted columnar defects using both RG and molecular dynamics methods. We find that the 
Splay glass transition temperature decreases for increased tilting of the columnar pins. The 
functional form for the transition temperature derived from an RG treatment is convincingly 
confirmed by our numerical simulations. We note that the experiments of Ref. [7] were only 
performed for small angle splay, vd ~ 0.176, and thus it would be a useful check of the 
theory to repeat these experiments with larger degree of splay. A full finite size scaling 
analysis is needed to determine the critical properties of the Splay glass transition. 
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Devereaux et al, Fig. 1 



FIG. 1. Splay glass transition temperature as a function of the width V£> of the tilting angle 
distribution and the elastic anisotropy parameter a = eia/uo, for various values of A x /K x . The 
top, middle, bottom surfaces are for A x /K x = 0.2,2, and 20, respectively. 
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Devereaux et al, Fig. 2 

FIG. 2. Average velocity of flux line motion versus applied force for T=0.01 (circles), 0.1 
(squares), 0.26 (diamonds), 1.26 (up triangles), 1.81 (pluses), and 1.91 (down triangles). The 
solid line represents Ohmic behavior (slope of 1) and serves as a guide to the eye. The critical 
temperature is estimated to be 1.87. Runs were performed for Vd = 2 and averaged over 10 
disorder configurations. Here the statistical fluctuations are smaller than the size of the points and 
a = 4 is kept fixed. 
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Devereaux et al, Fig. 3 

FIG. 3. Results for the Splay glass phase transition temperature derived from molecular dy- 
namics simulations of Eq. (11). The circles represent runs for the tilting angle distribution widths 
vd = tan($), while keeping a = 4 fixed. 
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